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ABSTRAGT 

We report the results of a search for large velocity width, low-intensity line wings — a commonly used 
signature of molecular outflows — in four low redshift (ultra)luminous infrared galaxies (U/LIRGs) 
that appear to be dominated by star formation. The targets were drawn from a sample of fourteen 
such galaxies presented in Ghung et al. (2011), who showed the stacked CO spectrum of the sample to 
exhibit 1000 km s“^-wide line wings. We obtained sensitive, wide bandwidth imaging of our targets 
using the IRAM Plateau de Bure Interferometer. We detect each target at very high significance but 
do not find the claimed line wings in these four targets. Instead, we constrain the flux in the line wings 
to be only a few percent. Casting our results as mass outflow rates following Cicone et al. (2014) 
we show them to be consistent with a picture in which very high mass loading factors preferentially 
occur in systems with high AGN contributions to their bolometric luminosity. We identify one of our 
targets, IRAS05083 (VII Zw 31), as a candidate molecular outflow. 


1. INTRODUCTION 


Galactic winds (|Veilleux et al.l 1200511 can carry gas, 
metals, and dust from a galaxy to its halo, deplete the 
future fuel for star formation, and inject energy and 
momentum into the inter- and circum-galactic medium. 
Winds seen in neutral gas via absorption lin es are a near- 
ubiqu itous feature of star-forming galaxies (|Rupke et al.l 
1200511 . In galaxies with strong active galactic nuclei 
(AGN), there i s also evidence for widespread molecu¬ 
lar winds (e.g .. iFeruglio et al.l 1201 01 lA atalo et al.l [20111 : 


IVeilleux et alJl^l3t iCicone et al.l 12014 . 12015^ with sug¬ 

gestions of a correlation between AGN strength an d 
the amount of outflowing material (iCicone et al.l[20T^ . 
There is also evidence for molecular outflows, or at least 
material blown out of the disk, in galaxies dominated 
by star formation (e.g. dWalter et al.ll200^ iBolatto et al.l 


120131 : iSakamoto et ahlfeoi^ . However, the number of 
very powerful star formation-driven molecular winds re¬ 
mains substantially lower than the set of molecular winds 
believed to be driven by AGN. 

An excellent prospect to build a large sample of star 
for mation-driven mole cular winds comes from a result 
by iChung et al.l (j201lL hereafter Gil) based on obser- 
vations of ul t ralurn in ous infrared galaxies (ULIRGs) by 
IChung et al.l (|2009ll . IChung et al.l (|2009ll observed CO 
emission from 29 ULIRGs using the wide-bandwidth 
Redshift Search Receiver instrument, at the time 
mounted on the 14-m Five Golleges Radio Astronomy 
Observatory. Cll stacked this sample into a single 
high signal-to-noise spectrum and noted that this spec¬ 
trum showed very wide (~ 1000 km s“^), line wings 
that acconnted for, on average, ~ 25% of the total 
CO emission from the galaxy. This is exactly the sort 
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of signal expected from a powerful outflow of molecu¬ 
lar gas from a moderately inclined galaxy. It resem¬ 
bles t he signatures seen in individu al sources like NGC 
1266 (lAlatalo et ^1201111 . Mr k 231 (IFeruglio et alJl20l^ 
iMaiolino et al.ll201^ . and the ICicone et al.l (|2014[) sam¬ 
ple of galaxies. The result of Cll is especially inter¬ 
esting because they showed that the stacked signal ap¬ 
peared to come preferentially from galaxies whose spec¬ 
troscopic classification was either “Starburst” or “H II 
region” rather than Seyfert or LINER. Thus, the Cll 
sample offers one of the best prospects for building a 
sample of powerful molecular winds driven primarily by 
star formation. 


2. OBSERVATIONS 

Based on the claim by Cll, during the summer of 
2014 we used the IRAM Plateau de Bure Interferome- 
tei0 to observe four members of the Cll Starburst/H II 
subsample. We targeted IRAS05083-I-7936 (105083), 
IRAS10035-b4852 (1100035), IRAS17132-b5313 Jl7132/ 
and IRAS17208-0014 (117208), which Cll classified as 
bright (105083,117208) or intermediate (1100035,117132) 
in CO flux. Table [T] lists these targets and their ba¬ 
sic physical parameters. We observed each target in the 
PdBTs most compact configuration for one track (ss 4 
hours on source), targeting the ^^CO (1-0) line using the 
WideX correlator in its maximum bandwidth configura¬ 
tion. The resulting data have synthesized beams of 5.0 x 
4.7" (105083), 6.4X4.4" (110035), 4.7x3.8" (117132), and 
5.6 X4.4" (117208) and noise of 0.9 (105083), 1.4 (110035), 
0.9 (117132), and 1.5 (117132) mJy beam“^. The data 
were imaged and deconvolved with a channel width of 
25 km s“^ following standard procedures in GILDAS. Be¬ 
fore any other analysis, we subtracted the continuum in 
the image plane using line free channels to either side of 
the CO line and not including the region of the CN lines, 
which are the only other bright lines in our bandpass. 
The peak 3mm continuum flux density in our targets is 

® This work is based on observations carried out with the 
IRAM Plateau de Bure Interferometer. IRAM is supported by 
INSU/GNRS (France), MPG (Germany), and IGN (Spain) 
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Table 1 

Observed targets 


Galaxy 

loglO tjTIR 
(logic Lq) 

/agn 

SFR 
M© yr-l 

Geo) 

IRAS05083-I-7936 

11.9 

< O.l'^ 

90 

0.054 

IRAS10035+4852 

12.0 

< 0.33 

~ 100 

0.065 

IRAS17132-I-5313 

11.9 

~ 0.17 

60 

0.051 

IRAS17208-0014 

12.4 

< 0.1 

~ 250 

0.043 


Note. — IR luminosit i es fro m lbanders et a 0 except 

110035 from Hwang et al.l J200711. AGN fraction estimates from 
IFarrah et al. ~ l|2003l 'l. see also IKoss et al.l l|2013l’) for a hard A- 
Ray upper limit on 117208. SFR translated from IR following 
ICicone et al.l J2014I) for consistency. 

^ Adopted based on high 6.2/xm equi v alent w idth of 0.64 (highly 
star forming) from IStierwalt et al.l For comparison, 

117208 has 0.31. 

~ 0.8 ± 0.04 mJy beam“^ (105083), < 0.36 mJy bearn”^ 
(5cr, 110035), and 0.63±0.03 mJy beam-^ (17132), 6.3± 
0.17 mJy beam“^ (117208). For 105083 and 11713, these 
are consistent with a synchrotron spectral in dex and the 
fluxe s of a few mJy measured at 33 GHz (iLerov et al.l 
1201 If) . For 117208, the 6 mJy flux is slig htly lower than 
the ft! 9 ± 1 mJy measured at 85 GHz by llmanishi et al.l 
()2006f) , but consistent within likely calibration uncertain¬ 
ties, especially if the emission has a synchrotron spectral 
index. 

A few residual sidelobes remain present in channels 
containing the main CO line for the two brightest tar¬ 
gets (105083 and 117208; see the peak intensity maps in 
Figure [T]). These could be removed with more aggressive 
cleaning but do not matter to the results of this paper. 
The bright line is already almost entirely cleaned and 
we search for line wings in channels away from the main 
line, where these sidelobes are not present. The dynamic 
range in the images is still very high, 350 in the highest 
case (105083). 

In this most compact configuration, the PdBl is ex¬ 
pected to resolve out sources larger than about 15" in 
an individual channel. This is much larger than the 
visible extent of any of our targets and the correspond¬ 
ing physical distance (« 15 kpc) is larger than the ex- 
tent of any known m olecular outflows at a fixed velocity 
(|Gicone et al.ll201^ . In our observations we find 96, 93, 
82, and 130% of the flux in the Cll FCRAO spectrum 
and the spectral shapes match well for the bright part 
of the line (see below). Given this good match to the 
single dish data and the large distance to our sources (so 
that even large physical scales are compact on the sky; 
15" Ri 15 kpc), we do not expect spatial filtering by the 
interferometer to be an important consideration for this 
work. 

T o calcul a te C O luminosities and sizes, we fol¬ 
low iWrightI (|2006ll using standard parameters {Hq = 
70 km s“^ Mpc“^, Hm = 0.3, Ha = 0.7; consistent with 
the references used for luminosity). For these parameters 
and a typical z = 0.05 the luminosity distance is 222 Mpc 
and angular size distance is 202 Mpc, so that 1" ~ 1 kpc. 

3. RESULTS 

We detected bright ^^CO emission from all of our 
sources. For each target. Figure [I] shows the inte¬ 
grated intensity (“moment 0”) maps (far left), intensity- 
weighted velocity (“moment 1”) maps (middle left), and 


maps of the peak intensity both along each line of sight 
as a function of sky position (middle right) and collapsed 
along the right ascensior0 axis (far right). The peak in¬ 
tensity maps saturate at 20 mJy beam“^ in order to show 
extended structure. In the first three panels, boxes show 
the region of interest encompassing bright emission for 
each target. We use this region to construct an inte¬ 
grated spectrum for each target. 

In two sources, 1100035 and 117132, we resolve the 
source into multiple components. In 110035 the two 
sources each show a velocity gradient that is likely rota¬ 
tion and so appear to be individual galaxies. In 117132, 
the second component is kn own to be a second , fainter 
galaxy from optical imaging (lArmus et al.lll990ll . 117208 
also shows some faint extension to the northwest, but 
this appears approximately consistent with an extension 
of the velocity field of the rotating disk. This system does 
have two components fsee lGarcia-Burillo et al.ir2015fl but 
they are not resolved in our observations. Therefore we 
treat this target as having only a single component. 

Table m reports the integrated properties of the CO 
line in each component. We integrate fluxes over the 
range ±1000 km s“^ around the nominal redshift, report 
the width of CO emission at 20% of its peak value (indi¬ 
cated by the dashed lines in Figure [5]), and convert the 
line flux to a luminosity using the distances described 
above. The line widths are reasonable for massive disks, 
200 -600 km s"^, and c omparable to other ULIRGS (e.g., 
see lChung et al.l 1200911 . Gradients, which are most nat¬ 
urally interpreted as signatures of rotation, are clearly 
visible in the velocity field maps, though with 5 kpc 
resolution, we only marginally resolve each target. The 
luminosities are typical of ULlRGs and consistent with 
previous measurements (see above). 

3.1. Search for Faint Line Wings 

The goal of our observations was to search for the pres¬ 
ence of faint line wings. Such wings are not immediately 
apparent from the images and position-velocity diagrams 
in Figure [H To search more rigorously, we collapse each 
cube into a single spectrum, integrating over the region of 
interest indicated in Figure [TJ We show these spectra in 
Figure [2] and report their properties in Table [2] For each 
spectrum, we plot two views: a version showing the full 
line at 25 km resolution and a version at 50 km s“^ 
resolution zoomed in to a y-axis range sa 10 times the 
rms noise in the spectrum. The velocity scale in each 
spectrum is relative to the nominal redshift of the source 
(the slight offset of IRAS05083 from zero indicates that 
our adopted redshift for the observations was slight off 
from the true redshift of the galaxy). 

Each source shows broad, bright lines but our obser¬ 
vations do not suggest the presence of the wide, faint 
component seen by Cll. In their stacked spectrum, the 
line wings have average magnitude r: 10% of the peak 
of the spectrum and contribute R 25% of the total CO 
flux. In each of our individual targets, a wide compo¬ 
nent 10% at 10% of the peak would be several sigma in 
each 25 km s“^ channel and extend across ^ 20 chan¬ 
nels. Such a signal is clearly not present in our data. If 

® We choose to measure the position-velocity diagram along the 
declination axis because our targets (conveniently) show velocity 
gradients mostly lined up in this direction. 
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Figure 1. CO observations of our targets, {far left) Integrated intensity maps. Gray contours mark regions with SNR > 3 (solid) and 
< —3 (dashed). Boxes in the first three panels show regions of interest used to construct the spectra used below, {middle left) Intensity 
weighted velocity field constructed from emission with SNR> 5. Velocity gradients, possibly indicative of rotation, are visible in all four 
targets, {middle right) Peak intensity maps for each target, saturated at 20 mjy beam“^ with contours at 3 (t and increasing by factors 
of 2. {far right) Peak intensity measured along the declination axis (which aligns roughly with the velocity gradient) within the region of 
interest. A circle indicates the beam extent in the left two images. Coordinates are plotted in each panel as arcsecond offsets from the 
pointing center; typically 1" ^ 1 kpc. Offset (0, 0) corresponds to the pointing center. 110035 and 117132 are each resolved into two discrete 
components. The offset structures around 105083 and 117208 are residual sidelobes in the main bright part of the line and not relevant to 
the current analysis. 


the line wings claimed by Cll are real, then they may be 
present in the other 10 targets that contribute to their 
stack of star forming galaxies. However, we can confi¬ 
dently reject the signal in four of the brightest members 
of this sample. 

Columns 6 and 7 of Table [2] quantify the flux in faint 
line wings in our targets. We report the integrated flux 
in the spectrum that is within ±1500 km s“^ of the nom¬ 
inal redshift of the source but outside the velocity range 
where the main line is > 5% of its peak amplitude. This 
region appears as the gray area in Figure [2j We report 
the flux in this velocity range (column 6) and express this 
as fraction of the total line flux (column 7). Only a few 
percent of the flux in our targets emerges from faint line 
wings, in stark contrast to the results found for a larger 
sample by Cll. In Section |4l we discuss this flux in line 


wings in terms of H 2 mass and mass outflow rate. 

3.2. Stacking Comparison 

We explore the contrast with Cll more directly in Fig- 
ure[3l The right panel shows a direct comparison between 
the spectra from Cll and the PdBI. For the bright part 
of the line in individual sources, the two telescopes agree 
quite well in both amplitude and line shape. In the left 
panel, we plot the stacked PdBI spectra for our four tar¬ 
gets. We follow a procedure similar to Cll: we first 
smooth the spectrum to a resolution of 125 km s“^. Then 
we identify the peak of the spectrum. We renormalize the 
peak amplitude of spectrum to unity and recenter the ve¬ 
locity axis so that the peak lies at zero velocity. We grid 
the spectra onto a grid centered at zero with channels 
125 km s“^ wide. Then we coadd all four spectra to 
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Table 2 

CO Properties From Integrated Spectra 


Galaxy 

Noise^ 

(mjy) 

,,,COb 

"’20 

(km s ^) 

Leo 

(Jy km s“^) 

Leo 

(10® K km pc-^ 

pwingsc 

^CO 

) (Jykms 1) 

rwingsd 

JCO 

IRAS05083 

2.0 

250 

100.2 ±0.4 

13.5 ±0.06 

1.1 ±0.5 

0.011 ±0.005 

IRAS10035 

2.4 

400 

38.4 ±0.5 

7.5 ±0.1 

0.9 ±0.6 

< 0.045 

IRAS17132 

1.9 

550 

54.3 ± 0.4 

6.6 ± 0.05 

0.1 ±0.5 

< 0.026 

IRAS 17208 

4.3 

550 

140.5 ± 1.0 

12.2 ±0.08 

3.2 ±1.0 

0.026 ± 0.007 

Note. — ! 

Statistical uncertainties 

are reported 

in the table, but the 

main uncertainty 

on 1020, ^CO, 


and Leo is systematic, not statistical. We adopt 15% plus uncertainty in the Hubble flow for Fqo 
and Loo km s~^ for W 20 . 

^ Noise in the integrated spectrum using 25 km s“^ channels. 

^ Velocity width at 20% of the peak value of the CO line. 

^ Flux in a broad, low component obtained by integrating the spectrum where |u| < 1500 km~^ and 
the flux is less than 5% of the maximum value. This is the gray area in Figure [2] 

^ The flux of the broad component expressed as a fraction of the total flux. For a Gaussian profile, one 
expects /wings ~ 0.014. Limits are 3a for IRAS10035 and IRAS17132. 
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Figure 2. Spectra of our sources integrated over the regions of interest shown in Figure^ Each spectrum is shown twice: first showing 
the full CO line {left, using 25 km s“^ channels) and then zooming in to focus on faint emission near the level of the noise {right, using 
50 km s“^ channels). The gray region shows the part of the spectrum used to measure the flux in line wings. Dashed vertical lines show 
W 20 i the width of the profile at 20% of its peak value. Table [2] reports integrated quantities for each spectrum. 
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Figure 2. (continued). 
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Figure 3. (left) Stacked spectra, normalized to the peak of each spectrum from our PdBI observations (blue) over stacked FCRAO 
spectra of the same four targets from Cll (gray) at 150 km s“^ spectral resolution. As discussed in Cll, the stacked FCRAO spectrum 
exhibits broad line wings at 10% of the peak flux density. This is evident even in the FCRAO stack of only our four targets. However, 
the PdBI spectra of the same galaxies, which have an order of magnitude better sensitivity, do not detect these wings or corresponding 
features in individual spectra (blue) despite a very good match of the observed (main) lines otherwise. The difference between the stacked 
spectra from FCRAO and the PdBI is shown in the bottom panel. This highlights the presence of a faint, broad component at ~ 10% of 
the peak intensity in the FCRAO spectra but not the PdBI spectra, (right) Comparisons of individual spectra. The PdBI observations 
match the FCRAO results for the bright part of the line, but the faint line wings implied by the FCRAO stack are absent. 


produce the averages. We carry out the same procedure 
for both the PdBI spectra and the FCRAO spectra. We 
plot both stacks in the left panel of Figure[3]and show the 
difference of the two stacks in the bottom panel. Note 
that the asymmetric frequency coverage around the line 
in the FCRAO spectrum leads to variable noise across 
the velocity range. 

The averaged FCRAO spectrum for our four targets 
shows a broad 1000 km s“^) component with intensity 
about 10% of the peak value. This closely resembles the 
stacked signal presented for 14 galaxies by Cll. That is, 
the stacked FCRAO spectrum for only our four targets 
does appear consistent with the claim of Cll for their 
entire star-forming sample. 

By contrast, the stack of the PdBI spectra do not show 
such a feature, so that the residual FCRAO minus PdBI 
plot in the bottom left panel of Figure [3] shows a signal 
that resembles the Cll claim near zero velocity. Rig¬ 
orously, we can only say tha t our observations disagree 
with the lChung et al.l (1200911 spectra regarding faint CO 
wings for these four targets. We cannot rule out the 
presence of very bright, broad line wings in the other 
10 targets. However given the disagreement for our four 
present targets we believe the current observations cast 
some doubt on the broader Cll claim. One possible ex¬ 
planation could be that baseline issues with the single 
dish telescope (which affect interferometers less) created 
a spurious signal in the single dish stack and that bad 
luck lead this signal to not appear in Oil’s carefully con¬ 
structed control samples. 


3.3. Faint Broad Wings 

Though we do not find a high-flux extended compo¬ 
nent, in both 105083 and 117208 we do detect line wings 
containing at most a few percent of the flux of the main 
line at marginal significance. These are visible in Fig¬ 
ure [5] and in the integrated line wing fluxes presented 
in Tabled! Using higher spatial resolution observations, 
IRAS 172 08 has recently been demo nstrated to have CO 
outflow (iGarcia-Burillo et al.l 120151 1 and the magnitude 
and velocity of the outflowing CO emission in their data 
(~ —500—750 km s“^ and a few mJy) is consistent with 
the faint emission in our spectrum). To our knowledge, 
105083 (VII Zw 31), does not have a previous measure¬ 
ment of faint line wings. Higher resolution observations 
will be needed to localize the high velocity emission rel¬ 
ative to the rotation curve and test whether the wings 
indeed reflect a CO outflow. 

4. DISCUSSION 

We fail to detect the line wings claimed by Cll, which 
account for rs 25% of the total line flux. Instead, we 
constrain the flux in faint line wings in our targets to be 
at most a few percent. How do these results, and the 
claim of Cll, compare with the r esults of the molecu - 
lar outflow population analysis in iCicone et al.l (l2014fl ? 
Qualitatively, that paper argued that high mass loading 
factors (defined as Mqf/SFR, see below) correlate to 
high contributions of AGN to the bolometric luminosity, 
/agn = UAGN/Uboi- In order to compare the various re¬ 
sults, we adopt the formalism and data compilation from 
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Figure 4. Molecular outflows, star formation, and AGN activity following the ICicone et all <120141^ formalism for Mqf (Equation [T)l. 
(left) Mass outflow rate as a function of star formation rate. A solid line marks equality, a mass loading facto r of unity. Dashed lines show 
ratios of 10, 100, and 1000. Our new points appear as red circles with black outlines. The ICicone et al.ril2014|) extended sample appears as 
gray circles marked in blue or red depending on their spectroscopic classification as “Hii” or “AGN.” Our estimate of the average properties 
implied by the claimed Cll stack appears as a large red square, (right) The mass loading factor, i.e., the mass outflow rate per unit 
SFR, as a function of AGN contribution to the total bolometric luminosity. Points and lines are the same as in the left panel. All of our 
measurements suggest mass loading factors < 10. None of the “Hii” targets in ICicone et al.l II2014I) show a high mass loading factor either. 
By contrast, reasonable assum ptions regarding the Cll stack would imply an average mass loading factor of ^ 30, something found only 
in AGN-dominated systems in ICicone et al.l ()2014I V The observations present here argue against such a high average mass loading factor 
in star-formation dominated galaxies. 


ICicone et al.l (l2014li . They assume a continuous outflow 
and calculate the mass outflow rate via 


,y q -^OF 

Mqf = 3u—— 

-KOF 


( 1 ) 


with V the mean velocity of outflowing material, Mqf 
the mass associated with the outfl ow, and Rg-p the ex - 
tent of the outflow. As discussed bv ICicone et al.l (I2014D . 
this geometry is likely to yield a high estimate for Mqf 
and we suggest that the reader bear this potential bias in 
mind. The advantage of this approach is that by adopt¬ 
ing it, we can make quantitative co mparative stat e ments 
between our measurements and the ICicone et al.l (|2014ll 
“extended sample.” 

We do not have an outflow extent for any of our tar¬ 
gets or the Cll sample and so ado pt the median Rg-p = 
600 pc from the extended sample of ICicone et al.l (1201411 . 
This is similar to the i?oF = 500 pc that they use for 
upper limits. For our targets, we adopt v ss 400 km s“^, 
approximately the Itr velocity for a 1000 km s“^ line 
width. The same rms velocity is a good description for 
the faint wings in 105083. For the Cll stack, we adopt 
{v) ~ 500 km s“^, half of the full extent of their lines, 
an a reasonable average velocity. 

We convert from lum inosity to Mqf using a 
iDownes fc SolomonI (|1998ll conversion factor of ago = 
0.8 Mq pc“^ (K km s“^)“^, again largely adopted for 
consistency. This implies Mqf < 1“3 x 10® M0 for 
our targets. We treat the measured fluxes for 105083 
and 117208 as upper limits on Mqf, reflecting the un¬ 
certainty in identifying them as more than outflow can¬ 
didates based on lower resolution spectra. Using this 
approach our targets have Mqf ^ 1“3 x 10® M0, which 


is similar to t he ty pical molecular mass for one of the 
ICicone et al.l (|2014ll outflows. For the Cll targets we 
take the average luminosity of the starburst subsample 
and the stacked result of 25% of the flux in the broad 
line. This yields an average Mqf ~ 1-5 x 10® M0. 

We calculate star formation rates from IR lumi nos- 
ity using the same formula as ICicone et al.l (I2014D for 
consistency. For the Cll stack, we adopt (SFR) ^ 
100 M0 yr“^, corresponding to a typical total IR lumi¬ 
nosity of 10^^ L0. Table [T] gives AGN fractions or limits 
for our targets, though we note that these are highly un¬ 
certain. We will discuss Cll’s starburst sample in terms 
of an upper limit of /agn < 0.5. 

Following this approach, our sources have upper lim¬ 
its < 150-400 M0 yr“^, which is < 1-4 times the star 
formation rate. For the Cll sample, the implied aver¬ 
age ou tflow rate is ~ 3000 M0 yr“^ in the ICicone et al.l 
(120141) formalism, which is « 30 times the SFR of an 
average individual target. 

As emphasized, the absolute values of these numbers 
is very approximate. Comparative measurements offer 
more insight and to this e nd we plot the resul ts of these 
calculations along with the ICicone et al.l (I2014D extended 
sample in Figure S) The left panel compares the mass 
outflow rate to star formation rate, so that diagonal 
lines correspond to mass loading factors (Mqf /SFR). 
The right panel compares the mass loading factor to the 
fraction of luminosity in the system supplied by an AGN. 
In both panels, blue points show galaxies spectroscopi¬ 
cally classified as an AGN and red points show starburst 
(“Hii”) systems. The three samples appear as different 
symbols. 

Considering only our points and ICicone et al.l (I2014D 
(i.e., leaving aside Cl I), the two plots reinforce the con- 
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elusion that mass loading factors above ~ 10 are only 
achieved in systems with a strong AGN contribution. 
The exact mass loading factors, especially for the star 
forming systems, may be slightly higher than the val¬ 
ues shown — for consistency all points here assume the 
same aco = 0.8 Mq pc“^ (K km s“^)“^. The geome¬ 
try may also vary, with poor resolution and inclination 
effects potentially hiding some outflows. Nonetheless, 
the data suggest that star formation appears to produce 
mass loading factors of up to a few, rather than ^ 10-100. 
Though we note the lack of a large set of high star forma¬ 
tion rate, SFR-dominated targets. By contrast, the Mqf 
suggested on average by Cll imply a mass loading factor 
of ^ 30 for their stack of starburst galaxie s. This would 
be as high as the strongest AGN in the ICicone et al.l 
(|2014ll extended sample. 

5. SUMMARY 

The calculations leading to Figure 0] require a large 
number of assumptions, but the figures still summarize 
the results of this paper. We follow up a claim of ex¬ 
tended, high-flux line wings in starburst galaxies by Gil. 
If present, these line wings would imply high mass load¬ 
ing factors in star-formation dominated system, as high 
as those seen in extreme AGN-dominated systems. We 
fail to find such features in 4 of the 14 galaxies pre¬ 
sented by Gil. Instead we limit the flux in such an 
extended spectral component to less than a few per¬ 
cent. We show that the stacked Cll signal is present 
in their data for our targets but not in the ^ 10 times 
more sensitive PdBI observations. Two of our targets do 
show marginally significant flux in an extended line wing: 
117208 is kno wn to have an outflow (lik ely driven by 
a hidden AGN iGarcia-Burillo et "ni 120151) . We identify 
the other, 105083 (VII Zw 31), as a candid ate molecular 
outflo w. Combining these results with the ICicone et al.l 
(|2014D extended sample, our observations of high SFR, 
star-formation dominated systems reinforce the idea that 
high mass loading factors tend to arise in systems where 
AGN contribute a large fraction of the luminosity. In 
addition to follow-up on our candidate outflow, sensitive 
observations of a large set of star-formation dominated 
high SFR systems will be useful to further test this con¬ 
clusion. 
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Table 3 

CN 1-0 Observations 


Galaxy 


113.5 GHz 
(Jy km s~^) 

113.2 GHz 
(Jy km s“^) 

CN113.5 

CN113.2 

CN113.5 

COl-0 

105083 

110035 


5.1 ±0.4 

3.6 ±0.4 

1.4 

20 

... comp. 

A 

0.63 ±0.24 

0.53 ± 0.24 


^ 30 

. .. comp. 
117132 

B 

0.87 ±0.32 

0.47 ± 0.32 


^ 25 

... comp. 

A 

1.6 ±0.3 

0.86 ± 0.33 

1.9 

30 

. .. comp. 

B 

0.58 ±0.13 

0.26 ± 0.13 


^ 10 

117208 


11.2 ±0.8 

5.0 ±0.8 

2.2 

13 


APPENDIX 

SERENDITIOUS CN MEASUREMENTS 

The wide bandwidth that allows us to search for CO line wings also allows us to detect the CN 1-0 J = 3/2- 
1/2,F = 5/2-3/2 and J = 1/2-1/2, F = 3/2-3/2 transitions in each target (rest frequencies 113.491 and 113.191 GHz). 
CN is a tracer of dense gas excite d in PDR regions; the line strength of the stronger 113.5 GHz line is typically within 
a factor of a few of HCN 1-0 (e.g.. lAalto et al.ll2002D . Table [3] reports the strength of the CN lines for each component 
and their ratios to one another and CO. The lines are at least tentatively detected for each spectrum. The ratio of 
the brighter line to the CO line varies between 10 and 30, consistent with previouswork. For the high signal-to-noise 
spectra and the two CN lines show ratios 1.4-2.2. The CN lines appear to have similar width to the CO, implying 
that the dense gas in these targets may be cospatial with the bulk of the molecular gas traced by CO. 










